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•Simple  analytic  approximate  expressions  have  been  derivt'd  t'L)r  1 
propagation  constants  in  the  earth-ionosphere  waveguide.  The 
expressions  are  an  extension  of  previous  w(jrk  to  a  more  gcuier a 
class  of  ionospheric-conductivity  height,  profiles.  'I’ln  express 
derived  here  allow  the  rapid  computation  of  EE!'  ))hase  sjua  ds, 
attenuation  rates,  and  excitation  factors  for  a  wide  t  itiae  f>f  1 
osphere  conditions  without  the  necessity  of  longtliy  ;ull-wa\’e 
computer  calculations.  The  results  can  be  aunlic'd  10  ;  he  t  a!>i^ 
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1 .  INTRODUCTION 

In  two  recent  papers  (Greifinger  and  Greifinger  [1978] 
and  Greifinger  and  Greifinger  [1979],  hereinafter  referred  to 
as  Paper  I  and  Paper  II) ,  the  authors  derived  simple  analytic 
approximate  expressions  for  TEM  eigenvalues  for  ELF  propaga¬ 
tion  in  the  earth-ionosphere  waveguide,  assuming  that  the 
ionospheric  conductivity  profile  could  be  approximated  in  a 
certain  convenient  way.  These  expressions  make  possible  a 
rapid  calculation  of  ELF  phase  velocities  and  attenuation 
rates  without  the  necessity  of  lengthy  full-wave  numerical 
calculations.  Comparisons  were  made,  for  a  number  of 
icnospheric-condt'ctivity  profiles,  between  approximate  and 
fall-wave  eigenvalues,  and  the  results  were  found  to  be  in 
excellent  agreement  in  all  cases.  More  recently,  Bannister 
[1979]  has  used  our  approximate  theory  and  the  Wait  [1964, 
1970]  VLF  exponential  ionospheric-conductivity  profile  to 
determine  propagation  constants  for  ELF  daytime  propagation 
in  the  earth- ionosphere  waveguide.  He  showed  that  the  result¬ 
ing  values  of  ELF  attenuation  rate  and  phase  velocity  are  in 
excellent  agreement  with  measured  data  from  various  sources  in 
L.ie  ircqucncy  range  from  5  to  200u  Hz. 

The  work  of  Bannister  clearly  demonstrates  the  validity 
and  ca Iculational  utility  of  tlic  approximate  theory.  It  also 
adds  creaence  to  the  belief  that  the  theory  should  be  equally 
applicable  to  PCA-disturbed  ionospheres,  for  which  comparable 
uat  a  acj  not  yet  exist,  and  to  nuclear-disturbed  ionospiieres . 

As  such.  It  is  [jotcntially  very  useful  for  making  rapid  and 
inexpensive  sensitivity  studios  of  ionospheric  pro;.)aqat ion 
under  disturbed  conditions.  It  '.herefore  seems  useful  to 
extend  the  approximate  theory  to  a  more  oeneral  class  of 
ionospheric-conductivity  iirofilcs  than  those  cons  Ldorc'd  in 
Papers  I  and  II.  The  conductivity  profiles  to  v;h  i  eii  tne 
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results  in  Papers  I  and  II  apply  will  appear  as  limiting  cases 
of  the  more  general  profile.  In  addition/  approximate 
analytic  expressions  will  be  derived,  within  the  framework  of 
the  theory,  for  the  excitation  factor,  which  is  an  important 
additional  measurable  quantity  in  ELF  propagation  experiments. 
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SUI4MARY  OF  PREVIOUS  WORK 


II . 


In  Paper  I,  the  authors  developed  the  theory  for  an 
ionosphere  that  was  sufficiently  disturbed  that  there  was  no 
significant  penetration  of  the  electromagnetic  field  to 
altitudes  where  anisotropy  due  to  the  Earth's  magnetic  field 
had  to  be  taken  into  account.  It  was  shown  that  the  propaga¬ 
tion  constant,  under  these  conditions,  was  determined  by  four 
parameters  —  two  frequency-dependent  altitudes  and  the  local 
conductivity  scale  height  at  each  of  those  altitudes.  The 
lower  altitude,  denoted  by  li^ ,  is  the  altitude  at  which  the 
conduction  current  becomes  equal  to  the  displacement  current. 
The  higher  altitude,  denoted  by  h^,  is  whore  the  rocii:'rocal 
of  the  local  wave  number  becomes  equal  to  the  local  scale 
height  of  tnc  refractive  index.  If  the  conductivity  [.rofilo 
in  the  neighborhood  of  h  and  h.  can  be  adccruatel'/  annu'oxin.at 
by  a  local  exponential,  tlie  expression  for  the  e  i  .jc  n\'. )  i  cl 

is 
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(  1  ) 


wlicre  ■  Q  and 
respectively 
velocity 


arc  the  local  scale  heights  at  iIq  and  h^, 

From  the  definition  of  S  ,  the  relativi'  i  ii.ist 
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and  the  horizontal  attenuation  rate  in  decibels  per  meejameter 


a  =  8. 68  X  10^  -  Im(S  ) 
c  o 


0.143  f 


hg  h^ 


wiierc  f  is  the  frequency  in  Hertz. 

In  Paper  II,  the  theory  v;as  extended  to  include  the  effects 
of  anisotropy.  For  the  sake  of  mathematical  simplicity,  it 
was  assumed  that  the  Earth's  field  was  vertical,  but  the 
results  actually  apply  over  a  rather  large  geographic  range. 
When  anisotropy  is  included,  the  parameters  which  determine 
tile  piropagation  constant  differ  for  daytime  and  nighttime 
ionospheric  conditions.  For  both  sets  of  conditions,  two 
parameters  which  enter  are  the  frequency-dependent  altitude 
Uq  at  wiiich  tiie  conduction  current  parallel  to  the  magnetic 
fieiu  becomes  equal  to  the  displacement  current,  and  the  local 
scale  iieight  of  the  parallel  conductivity  Under  daytime 

conditions,  two  additional  parrs  of  frequency-dependent 
altitudes  and  associated  scale  heiglits  aiipear  as  paramotoi's. 

One  altitude  is  that  at  wiiich  the  reciprocal  of  the  local  v;a\'i 
number  for  vertically  propagating  O  waves  becomes  equal  to 
tile  local  scale  height  of  the  refractive  index  for  such  waves. 
The  associated  scale  height  is  the  scale  height  of  tiiis  refrac¬ 
tive  index.  The  other  pair  of  parameters  are  the  corros j>ond i  iig 
quantities  for  vertically  propagating  X  (whistler)  wavxs.  7 
these  altitudes  arc  attained  in  a  region  of  the  ionosphi'ro 
wave  j  jjjl  l'-'p|/  where  Ojj  is  the  Hall  conductivity  anci  Liie 
E’edersen  conductivity,  the  two  pairs  of  parameters  becoiiK  iden¬ 
tical.  Tliis  is,  in  fact,  the  case  over  a  substantial  altitiuu' 
range.  Under  these  conditions,  tJie  app>roximatt'  exprt's  r,  i  cmi  1  (,ir 
the  I'EM  eigenvalue  is 


where  the  altitude  and  scale  height  are  referred  to  the 
Hall  conductivity.  This  expression  is  slightly  different  from 
Eq.  (1),  its  counterpart  for  the  isotropic  ionosphere.  How¬ 
ever,  since  <<  1  and  1/  the  value  of  is  the 

same  to  lowest  order  in  these  quantities. 

For  typical  nighttime  conditions,  a  sharp  reflecting 
E-region  bottom  may  be  encountered  before  the  reciprocal  of 
the  local  wave  number  becomes  equal  to  the  local  scale  height. 
Under  such  circumstances,  the  altitude  of  the  E-region  bottom 
replaces  h^^  as  a  parameter,  and  the  local  wavelength  on  the 
E-region  side  of  the  bottom  replaces  the  scale  height  as  a 
parameter.  Under  these  conditions,  the  approximate  TEM 
eigenvalue  is  given  by 


where  =  i./c  and  n^  is  the  index  of  refraction  tor  X  waves 
just  inside  the  E-region. 

The  approximate  daytime  expressions  were  deriv<_‘d  undvr  I 'r.v 
assumiition  that  the  conductivity  profile  in  tin.'  neiqhijoi  (imck; 
of  hj^  could  be  li'prcjx  imated  by  a  local  expcjnential  over  an 
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altitude  ranue  of  a  scale  heiLiiit  or  so.  Tiii..  apr^rox  iniate 
nightuinie  expression  was  derived  under  the  assumption  of  a 
sharp  transition  in  the  ionosp'hor  ic  conductivity  profile  do  J  ow 
wnich  the  local  wavclenqth  is  larqc  comp'arcd  witii  the  local 
scale  heiqht  and  above  which  it  is  small  comt.)arcd  with  tnu 
local  scale  height.  There  is  an  important  class  of  ionospiior i c- 
conauctivity  profiles  which  bridge  the  gap'  betv.'oen  tne  assumed 
daytime  conditions  and  nighttii-ie  conditions.  T'no'se  are  condi¬ 
tions  wh.cre  the  altitude  is  attained  in  a  ro'gion  of  rapid, 
but  not  discontinuous,  increase  of  conductivity  with  altitude. 
Under  such  circumstances,  the  conductivity  rofilc’  in  ‘ne*  nei- 
borhood  of  the  altitude  hj^  cannot  di-  accurat>.  ly  r.  presentee  by 
a  local  exponential.  In  the  followine  SL'.tii.n,  tiie  ap;  rox  i  ;v.a  ti. 
theory  of  Papers  I  and  II  will  be  extondeu  to  witii  sue;;  a 

situation. 


IT  I.  gi::;l!?alizl'd  co:a:)ucTiviTY  PkoriLi; 


In  this  section,  we  consider  a  more  ^  ionosphe’r i c- 

conduct iv i ty  profile  in  the  neiqhborhood  of  tiian  tiiostj 
treated  in  Papers  I  and  II.  As  in  the  (irovious  work,  is 
defined  as  the  altitude  at  which 


wiiero  is  the  conductivity  at  the  altitude  h^  anc 


1  do  1-1 


is  the  local  scale  heiqht  at  this  altitude.  A  more  qeneral 
conductivity  profile  in  the  vicinity  of  n^  can  bo  written  as 


( z-'n  )  / 

(1  +  s)  c 
1  (z- 

1  '  se 


whc-re  s  IS  an  additional  parameter.  Tin,-  sc.'ilc'  heioiit 
re  la  t  I'd  to  the  other  T>arai;\t.:tcr s  throuah 


d  (  1  +  s  1 


wiiLch  fcjJlows  from  i.es.  lo)  air,  (9). 


Tin.  ■  i  .1 1  e  nish  in  ■  :  i.  : .  (  'i ;  in.;  •  , 


:  .  1 !‘  iru  t  ■  I  s  as 


s 


(11) 


=  e 


where  Eq.  (11)  defines  the  altitude  .  The  profile  given  by 
Lq.  (9)  then  talcos  the  form 


(z-h  ) / • 

(1  +  s)  c _ 

1  (z-h  )/' 

1  +  e  ^ 


(12) 


This  profile  is  approximately  exponential  up  to  altitudes  a 
scale  heigiit  or  so  below  h  and  becomes  essentially  constant 
a  scale  height  or  so  above  hj,.  This  is  the  kind  of  transition 
'vvhich  typically  occurs  between  the  D-rogion  and  the  E-rogion, 
and  h^  may  thus  be  thought  of  as  the  lower  boundar}-  of  the 
E-region.  If  hg>>hj^,  the  profile  in  the  neighborhood  of 
is  a  simple  exponential  with  scale  height  '  ■  ,  which  was 

tne  daytime  profile  treated  in  Papers  I  and  II.  I  f  h^.  •  -  h^, 
the  E-region  is  encountered  before  the  altitude  is  attained. 
If,  in  addition,  the  scale  height  ^  is  sufficiently  small,  the 
transition  between  the  D-  and  E-rogions  is  very  sharp,  v,-l',ich 
was  the  nighttime  profile  treated  in  Paper  II.  The  generalized 
conductivity  profile  given  by  Eq.  (9)  thus  contains  the  '.iro- 
vious  profiles  as  limiting  cases. 
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DL'TERM1.\’AT10.\  OF  APPHOXIM/\TL  TPM  EIOFNVALLFh 


As  in  Paper  II,  we  assume  that  the  Earth's  magnetic  field 
is  vortical.  In  the  QL  approximation,  propagation  in  the 
vicinity  of  h^  is  governed  by  two  uncoupled  wave  equations. 

In  terms  of  quantities  Iq  defined  in  Paper  II,  these  can  bo 
written  as 


2^2  , 


0 


(13) 


where 


n"  =  ^  (i.-p  •  (H> 

The  upper  sign  corresponds  to  vertical  O  wave  propagation  and 
the  lower  sign  to  vortical  X  wave  (whistler)  propagation. 

Under  normal  daytime  conditions  at  the  upper  end  of  tiii,  CUE 
band  and  for  sufficiently  disturbed  daytime  or  nighttime 
conditions  at  all  ELF  frequencies,  tiie  altitude  h^  is  attained 
in  a  region  where  the  ionosphere  is  nearly  isotropic,  so  th-t 


O' 

at  low  [rec;uencies  and  norn' 1 1 
frequencies,  hj^  occurs  when 


'll  ^  0'  p  ''  -'o 


2  2 
n^  -  n_ 


Under  normal  daytime  conditions 
nighttime  conditions  at  all  ELF 


H 


so  that 


(IG) 


2 


n 


+ 


-H_!  _ 
O'  i 


It  will  be  assumed  that,  whichever  case  applies,  the  relevant 
conductivity  profile  (Uj^or  .-ji)  is  yiven  by  Eq .  (9).  Both 

cases  may  be  treated  within  the  same  framework  by  writiny 


n 


2 


i 


e 


(1  +  s) 
1  +  se 


(z-h^)/-2 


(z-h^)/-2 


(II) 


where 


+ 


(18) 


for  the  isotropic  case  and 


'  + 


0 


(19) 


for  the  anisotropic  case. 

The  wave  equation  may  be  solved  by  introduciny  the  variable 

( z-hj^) 

y  =  -se 

which,  with  Eqs.  (7)  and  (10),  transforms  Eq .  (13)  to 

,2 

+  L  _ 1 _ 

,,2  V  iy  4s  ( 1  +  s)  y  (y  -  1) 


-  0. 


l.'v) 


This  can  bo  recognized  as  the  generalized  hypergeonnctric 
differential  equation.  The  solutions  which  obey  tlie  radia¬ 
tion  condition  at  large  altitude  (y  ■■  -  ■  )  are  (Magnus, 
Oberhettinger ,  and  Soni  [1966]) 


1  +  21-  ,  ; 


where  F  is  the  generalized  hypergeometr ic  function,  and 


i  (s'  .  -  TT)  /2 
c _ 

2[s(l  +  s)]'^ 

These  solutions  are  valid  down  to  an  altitude  a  few  scale 
heights  below  h^  or  hj,,  whichever  is  smaller.  At  these 
altitudes  (y  -  0) ,  the  solutions  take  the  form  (Magnus, 
Oberhettinger,  and  Soni  [1966]) 


*  C  ,  [  z  -  G  ( r  ,  )  ] 


whore  C.  are  constants,  and 


0  ( ;• ) 


2'.  (1  +  ;  )  -  2:  (1) 


-  +  loa  s 


wlu/re  r  i  tlK'  logarithmic  do'i'ivative  of  the  gainiua  luiiotiei'. 


The  eigenvMlU)-’  is 


r.(  j  i  u  t  i  <  -nr. 
Ueund  aiul 


wi'.  ich  obey 
are  valitl 
SI  j  1  u  1 1  on^; 


determined  by  nialciiiiu]  I'.g.  (  .:  )  to 

tile  proper  bourui.u  y  cotKiitions  it  tin 
up  to  altitudes  a  t(.'W  scale  ia  i.;lit,--  li 
wei'i-  shov;n  in  Fasaa  1  1  to  Inv''  (  in  l-’t 


Tiles  I 


=  A  1(1  ^  C)  z  -  s; 


where  A  and  C  are  constants.  For  z  >>  h^,  the  integrand  in 
Eq.  (2S)  becomes  very  small,  and  little  error  is  made  by 
extending  the  range  of  integration  to  infinity.  Matching 
Eqs,  (23)  and  (25)  in  the  region  where  both  are  valid  then 
leads  to  the  condition 


G(e^)  +  G(ej 


which  is  the  generalized  expression  for  tiie  TF.r'  eigenvalue. 
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COMPARISON  WITH  PREVIOUS  RESULTS 


The  approximate  oxpressioi\s  derived  in  Papers  1  and  ii  can 
be  recovered  as  limiting  cases  of  Uq.  In  those  ;.a]jers, 

it  was  assumed  that  the  conductivity  profile  in  the  neigh¬ 
borhood  of  Hq  coulu  be  represented  by  a  local  expotientiai , 

i.  e.  , 


'  o 


(^-ho)/-o 


o 


(27) 


With  this  approximation,  the  integral  in  Eej .  (2i)  can  be 

evaluated  analytically,  giving 


1 


1  • 

2 


o' 


(2i) 


The  eigenvalue  can  thus  be  written  as 


f  (, 


+  '■ 


^  -  2 


(20 


whore 


/i;(  ^  )  f;( 


'FLe  limiting  forms  of  will  now  h<  ■  ueti  ir.'in.'d. 

a.  Daytime  ionospheric  c:onci  i  t  ions- - I  ti  tin  pi(\’i''U;.  wi  a  L  , 
U  was  assum'd  that,  lor  d^iytimi-  conditions,  lii'  ill  i  i  ic.. 


was  below  the  bottom  of  the  b-region  and  that  the  conductivit 
profile  in  the  neighborhood  of  h^^  could  be  approximated  !"■  an 
exponential.  As  discussed  in  Section  III,  this  is  the  limit 
hj^^.  ’  S'  -  1  of  the  more  general  conductivity  profile.  Ii. 

this  limit, 


i  ( 


Since  i  -  ,  |  1,  we  niay  ust-  thc'  asy.-^g.t  ( a  i  c  I'xpansion  of  the 

•  function  to  t^blain 


.  (1  +  .  )  5;  lo  j  - - - -  -  log  2  -  log  s.  {  b!) 

Ctimbined  with  l.g.  ( i  1 )  ,  this  gives 

C.  (  ,  )  ■  hj  -  ■  1  1  (  -  )  -  2{  .  -  log  2)1  (  :  ,  ) 

'Where  ,  is  buK'r's  e'onstant.  I'oi'  tlu'  isotrojoc  daytimi.' 
ionospiH'f’  ol  Paper  1  ,  =  ^  ,  so  that 

r(  _)  hj  *  S'  +  2( .  -  lea  2)j.  (  I) 

In  the  limit  s  1,  ,  ~  'I'iuis  ,  aeiart  I'liu  a  '.■n,a  I  1  eoric'C 

tion  term,  wiiich  was  neglected  in  I’apL-r  I,  tiu  e  i  :eii\-a  1  ic 
obtained  from  [ijs.  (/'')  and  (  o} )  ,  agrees  w  i  t  ii  i:<i.  (1). 

l  oi'  the'  anisotropiC  da'/timi,  ionospheit,  ol  I’a;  <  i  11,  ^ 

Til  us  , 
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whore  the  small  correction  term  has  been  neglected.  Since 
Si  this  gives  the  same  eigenvalue  as  Eq .  (4)  . 

b.  Nighttime  ionospheric  conditions--It  was  assumed  in 
Paper  II  that,  for  nighttime  conditions,  a  sharp  E-region 
bottom  was  encountered  before  the  altitude  h^  was  attained, 
and  that  the  conductivity  increased  slowly  above  the  E-region 
bottom.  As  discussed  in  Section  III,  this  is  the  limit 
h^  ■  ■  hj^,  s  >  '  1  of  the  more  general  profile.  To  examine  this 
limit,  it  is  convenient  to  write 
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which  follows  from  Eqs.  (7),  (9),  and  (10)  with 
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[n  the  limit  s  ■  •  1,  I,-',  ■  1  and  (1  +  .')  ~  ,  (1)  .  Furtliermori' , 


log  s  =  (h^  -  hj.,)  /  ■ 


(  3  H ) 


Tl  liien  follows  from  Eqs.  (tl)  anci  (  ^0)  that 
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where 
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To  lowest  order  in  e,  the  eigenvalue  obtained  from  Eqs. 
and  (39)  agrees  with  Eq.  (5). 


(40) 


(29) 


We  have  thus  demonstrated  that  the  results  of  Papers  I 
and  II  can  be  recovered  from  the  more  general  ionospheric 
conductivity  profile. 
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VI. 


APPROX  I  r-lATE  EXPRESSIONS  FOR  EXCITATION  FACTOR 


The  ELF  signal  strength  at  a  receiver  depends  not  only  on 
the  propagation  loss  between  transmitter  and  receiver,  but 
also  on  the  excitation  factor.  Predictive  calculations  of 
the  effects  of  ionosj'heric  disturbances  on  ELF  signal 
strengths  must  therefore  take  this  factor  into  account.  In 
Papers  I  and  II,  approximate  expressions  were  developed  only 
for  the  phase  speed  and  attenuation  rate.  In  this  section, 
the  theory  will  be  augmented  with  approximate  analytic  expres¬ 
sions  for  the  excitation  factor  as  well. 

For  the  calculation  of  the  excitation  factor,  it  is 
convenient  to  reformulate  the  theory  somewhat.  Following  a 
customary  procedure  in  ELF  calculations,  we  first  make  a 
Fourier  decomposition  in  the  transverse  plane  by  writing 
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where  the  hat  is  used  to  denote  the  Fourier  t runs  f  ont:.  NA 
next  introduce  a  rectangular  set  of  ^-dependent  unit  victois 
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and  resolve  all  Fourier  transforms  into  their  S-dependent 
components.  VVhen  rhis  is  carried  out,  Maxwell's  eejuations 
separate,  in  the  isotropic  lower  ionosphere,  into  two  si.-Ls  oi 
uncoupled  equations,  one  qoverninc  the  propagation  of  TL  fields 
and  the  other  of  TM  fields.  The  equations  for  the  TM  fiolus, 
which  include  tlie  TdM  mode,  arc 
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wnere  x  =  k  z  ajid 
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(4-5) 


If  we  introduce  as  a  variable  the  wave  impedance 


icB^ 


( 4  ^ . ) 


Lne  Til  field  equations  can  be  combined  into  tlio  single  e'Ciuaticm 


n‘^2^ 


wnif:h  is  <i  noni-ineai,  [  i  rst-ordt'r  d  i  f  f  t- rea^  1 1  a  1  (Hiuatieii  et  f  iu 
1-’  1  e  t:a  t  i  i  y  j  n  . 


2  0 


We  consider  a  small  horizontal  electric  dipole  of  strenqtn 
IL,  oi'iented  along  the  x-direction  (,=0),  and  located  at  a 
iioigiit  in  a  thin  uniform  "source  layer"  above  the  ground. 
Tne  solution  in  tne  source  layer  can  be  written  as 


(z) 


u  ILcS  C 
o  X 


2S (1+rR) 


iCz, 

e 


-iCz 


-iCz  , 


+  Re 


c 


1 C  z , 
re 


( 


where  z  is  the  greater  of  z  and  z  and  z  is  the  smaller  of 

s 

those  two  quantities,  and 
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The  quantities  K  and  r  are  reflection  coefficients  to  be 
determined  from  the  boundary  conditions  at  the  up5;er  and  lower 
boundaries  of  the  source  layer,  respectively.  Thest  quaoititi., 
can  be  expressed  in  terms  of  the  wave  impedances  at  tht:se 
boundaries  by  the  relations 
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where  we  liave  f^laced  the  source  layer  at  the  nic.una  i  ■'  0)  ei;a 

allowed  it  to  bt.'come  in  f  i  n  i  Les  inia  ]  ly  tinn. 


1  O' 


The  b''Uiiuaiy  conditions  recjui  rc’  ttie  o.- ni  t.  i  r.  u  i  t  y  ■; 
tan'iL-nt  la  i  electric:  and  maanet  i  c  {rclao.  ..it  all  ini'.  i‘, 


This  allows  us  to  <.‘valuate  ^j^(O)  from  t/ie  solutions  in  the 
sjrounu.  In  duriviny  the  approximate  oiyenvalues,  it  v;as 
assuiueu  that  tiie  qrounu  was  perfectly  conductiny.  Tne  finite 
conductivity  of  the  grounu  introduces  only  a  small  correction 
to  these  expressions.  However,  the  assumption  of  infinite 
ground  conductivity  cannot  be  made  in  deriving  tiiC  excitation 
factor  for  a  horizontal  electric  dipole  at  the  ground,  since 
tiieie  woulu  be  no  excitation  of  the  waveguide  at  all  under 
those  conditions.  The  solutions  in  the  ground  obeying  tne 
radiation  condition  have  the  form 
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where 
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It  then  follows  from  Eqs.  (lib)  and  (4j)  that 
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The  impedance  is  determined  from  the  solutions 

above  the  source.  The  w'av'-e  eciuation  foi'  the  impodarce  can  l.i 
written  as  the  intc’yral  equation 


=  Z^,(0) 
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tho  contr ibutioTi 


It  is  not  difficult  to  show  that,  for  7  ■  ■  h^, 
to  tho  inteiral  from  the  last  term  in  tne  squiue  l^rackots  is 
small.  I'nus , 


(  i  I  ) 


wnere  is  given  by  Ei-j.  (44  j  .  Furthermore,  if  x  ■-  ,  little 

error  is  luade  in  extending  the  range  of  integration  to  infin¬ 
ity.  The  wave  impediince  above  tho  source  thus  takes  t:u.  form 


Z  S  Z^(0)  +  +  k^iiQ 

(h|^  •  ■  z  ■  hj^)  (  ,  4 ) 

wiiere  n^  is  defined  by  tlie  equation 
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For  an  exponential  conductivity 
right-hand  side  of  Eci .  (28). 

From  the  solutions  (jiven  in 
near  the  im[jedance  somewliat 


profile,  hy  is  ejiven  by  Li.<- 

Papers  I  and  II  fui'  tin.'  reana 
below  h  takes  tho  f<.>rm 


=  i.-  .  t .  (  z  -  h  ,  ) 


where 


and  t  t  ■  ^  )  is  aiven  by  Eq.  (  ;o)  •  For  the  earticuiar  nrol'ile 

treated  in  j)rovious  work,  i  Fy^/  _)  reduces  to  Eqs.  (,;4)  ,  i  h,)  , 
or  (  ,  as  the  casu  may  be. 

By  requiring  tile  ieadrnq  terms  of  Eqs.  (44)  and  (S'  )  to 
aqree  in  the  rcqion  where  both  are  valid,  wc 

obtain 
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(0)  =  Fq  (S  hp  -  hj^)  .  ( '  ■  ) 
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It  should  be  noted  tiiat  at  this  point,  we  may  regard  S"  as  a?', 
undetermined  pararaotor,  rather  tiian  the  TEM  eigenvalue. 

With  Z.  (0)  and  Z,.(0)  now  determined,  wc  can  e.xpress  the 
reflection  coefficients  r  and  R  in  terms  of  the  [parameter  S  by 
means  of  Eqs.  (4'Ja)  and  (lyb).  This,  finally,  allows  us  to 
(..'xpress  and  in  terms  of  this  paraiiiotor.  The  cxpressio.n 
for  at  the  grounu  becomes 
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where  we  have  neglected  terms  of  in-Lative  oidi.r  n^,  ,  wliici; 

cjre  small  at  ELI'. 

It  now  remains  to  carry  out  tlie  Lndicat(.'U  i  n  f  ^  ;  i  a  t_  i  on  in 
Eq .  (4Ia)  .  Cariying  out  tne  integration  ovni  uniK  in  il-si'ace 

wc  obtain 
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In  the  remaininq  inteqration,  the  inteqral  can  be  extended  to 
-■  by  repLacincj  the  Bessel  function  by  an  appropriate  Hankel 
function.  The  intcyral  can  then  be  evaluated  by  contour 
integration  in  tiie  coniplex  S-plane.  At  large  distances  frojii 
the  source,  tiie  primary  contribution  is  f rom  t!ie  pole  of  flu,' 
integrand  at 
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( 0 ; ) 


uiiich  will  be  recognized  as  the  approximate  TL'M  eigenvalue 
obtained  previously.  This  may  be  considered  as  an  alternative 
derivation  of  that  expression.  From  the  residue  at  this  i-ole, 
we  obtain 


I  i  ,(.Xi.S  . 

'•  J/ 


e  e 


^n,.  I 

G  a  g 


( ' 


wl’.ere  ■  is  the  wava'leiiuth  in  1  law  space  and  wc  iia'.n  a.'c 
large  argument  expansion  ot  tile  Hankel  iunction. 

We  lielint^  tin  magnetic  excitation  tactm  Ijv  tin 

equation 
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which  differs  somewhat  from  otiier  definitioiu;  (e.q.,  I-iaiii; i s L(  r 
[1975])  in  separatirn-!  the  effect  of  tiie  ionosphere  fro;;;  ti'.af 
of  tlie  ground.  With  tJiis  defiiiition,  th<  excitation  facnji 
becomes 


,  _  ^-1/2  _  -3/4  -1/4 
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wnere  h^  is  given  by  bq.  (37)  .  Finally,  since  b  ^  ,  tin- 

analogous  excitation  factor  for  the  vertical  electric  fiela  is 


-  3/4  ,  -  1/4 
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It  was  [jointed  out  in  Pa[jer  II  tiiat.  tiu  :u  r  i  in  -.t.i!  r 
of  energy  flow  is  essentially  constant  up  to  an  altituue  , 
above  which  it  falls  off  very  rapidly  with  altitude.  It  is 
tnerefore  not  surprising  that  tjie  effective  waveguide  iieighf 
for  excitation  is  rougnly  h^,  rather  tnan  the  higher  reflect¬ 
ing  liL'ight 


VII  . 


COWCLUSI(,),\S 


Tiic  approximate  theory  developed  in  Papers  I  and  II  for 
the  determination  of  ELF  propagation  constants  iii  the  earth- 
ionospiiere  waveguide  lias  been  extended  to  a  rather  more  gen¬ 
eral  ionospheric-conauctivity  profile  tiian  those  considered 
previously.  The  results  of  Pajjers  i  ana  I]  arc  recovc  r^a  as 
limiting  cases  of  the  more  general  profile.  Sing-ie  analytic 
expressions  have  also  been  derived  for  the  excitation  factor, 
wliich  is  another  important  parameter  in  ELF  propLigation . 
AltJiuugn  approximate  excitation  factors  have  not  bc'en  comj  aieu 
witii  full-wave  calculations,  it  is  reasonauie  to  assume,  cm  t!i«. 
basis  of  tile  generally  excellent  agreement  between  afjproximati' 
and  full-wave  propagation  constants,  tliat  tiiey  represent  vi-ry 
good  appro.ximations .  The  simple  analytic  expressions  derived 
in  tills  work  allow  tne>  rapid  computation  of  piiase  speens, 
attenuation  rates,  and  excitation  factors  for  a  wide  range  of 
ionospheric  conditions  witiiout  the  necessity  of  lengthy  tull- 
wave  numerical  calculations .  Tlic  thc'ory  can  tiiercforo  be 
quite  useful  in  the  study  of  the  effects  of  a  variety  of 
ionospiieric  disturbances,  both  natural  and  artificial,  on 
LLF  communication  systems. 
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